Mathematical models for passive imaging 
I: general background. 

o . 

O . Yves Colin de Verdiere 

O ■ February 7, 2008 

O 

Abstract 

■ Passive imaging is a new technics which has been proved to be very 
efficient, for example in seismology: the correlation of the noisy fields be- 

' tween different points is strongly related to the Green function of the wave 

■ propagation. The aim of this paper is to provide a mathematical context 
for this approach and to show, in particular, how the methods of semi- 

^ classical analysis can be be used in order to find the asymptotic behaviour 

of the correlations. 
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Introduction 



Passive imaging is a way to solve inverse problems: it has been succesfull in 
seismology and acoustics [21 El EEB UB EE 1201 12H I2H] • The method is as follows: 
let us assume that we have a medium X (a smooth manifold) and a smooth, 
deterministic (no randomness in it) linear wave equation in X. We hope to recover 
(part of) the geometry of X from the wave propagation. We assume that there 
is somewhere in X a source of noise i(x,t) which is a stationary random field. 
This source generates, by the wave propagation, a field u(x, t) = (u a (x, t)) a=1 ... N 
which people do record on long time intervalls. We want to get some information 
on the propagation of waves from B to A in X from the correlation matrix 1 
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(equivalently 



C a ,b(t) = lim — / u(A,t) (g)u(B,t- t) dt 



Cf )B {r)= lim 1 f u a (A,t)uP(B,t-r)dt 
T ^°° 1 Jo 



Tnstitut Fourier, Unite mixte de recherch e CNRS-UJF 5582, BP 74, 38402-Saint Martin 
d'Heres Cedex (France); http://www-fourier.ujf-grenoble.fr/ycolver/ 



1 For every matrix (oy), we write {flij)* ■= {a>ji) 
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which can be computed numerically from the fields recorded at A and B. It 
turns out that Ca,b{t) is closely related to the deterministic Green's function 
G(A, B, t) of the wave equation in X. It means that one can hope to recover, 
using Fourier analysis, the propagation speeds of waves between A and B as a 
function of the frequency, or, in other words, the so-called dispersion relation. 

If the wave dynamics is time reversal symmetric, the correlation admits also 
a symmetry by change of r into — r; this observation has been used for clock 
synchronization in the ocean, see [T7j . 

The goal of this paper is to give precise formulae for Ca,b(j) in the high 
frequency limit assuming a rapide decay of correlations of the source f . More 
precisely, we have 2 small parameters, one of them entering into the correlation 
distance of the source noise, the other one in the high frequency propagation. 
The fact that both are of the same order of magnitude is crucial for the method. 

Let us also mention on the technical side that, rather than using mode de- 
compositions, we prefer to work directly with the dynamics; in other words, we 
need really a time dependent rather than a stationary approach. Mode decompo- 
sitions are often usefull, but they are of no much help for general operators with 
no particular symmetry. 

For clarity, we will first discuss the non-physical case of a first order wave 
equation like the Schrodinger equation, then the case of a more usual wave equa- 
tions (acoustics, elasticity). 

The main result expresses, for r > 0, Ca.b(t) as the Schwartz kernel of Q(t)oH 
where II is a suitable pseudo-differential operator (a \l/DO), whose principal 
symbol can be explicitely computed, and fi(r) is the (semi-)group of the (damped) 
wave propagation. It implies that we can recover the dispersion relation, i.e. the 
classical dynamics, from the knowledge of all two-points correlations. 

In order to make the paper readable by a large set of people, we have tried to 
make it self-contained by including sections on pseudo-differential operators and 
on random fields. 

In Section^ we start with a quite general setting and discuss a general formula 
for the correlation (Equation (J3J). 

Section El is devoted to exact formulae in case of an homogeneous white noise. 

In Section |31 we discuss the important property of time reversal symmetry 
which plays a prominent part in the applications and is also usefull as a numerical 
test. 

In Section 01 we introduce a large family of anisotropic random fields and 
show the relation between their power spectra and the Wigner measures. 

Section |3] contains the main result expressing the correlation in the case of a 
Schrodinger wave equation. 

Section El does the same in case of a wave equation. 

The short Section [7| is a problem section. 

Section IS] is a about a quite independent issue relative to correlations of scat- 
tered waves. 



2 



Finally, there is a very short introduction to pseudo-differential operators in 
Section 

Acknowledgment: I would like to thank Michel Campillo and his colleagues 
from L GIT for discussions and collaborations; thank you also to Raoul Robert 
for the starting ideas about constructing random fields. 

1 A general formula for the correlation 
1.1 The model 

We will first consider the following damped wave equation: 

• X is a smooth manifold of dimension d with a smooth measure \dx\ 

• u(x, t), x G X, t G M is the field (scalar or vector valued) with values in 

(or R N ). 

• The linear operator H is the Hamiltonian, acting on L 2 (X, <C N ). It satisfies 
some attenuation property: if we define the semi-group Q(t) = exp(—tH),t > 
0, there exists k > 0, so that we have the estimate ||fi(t)|| = 0(e~ fc< ). 

• The source f is a stationary ergodic random field onlxl with values in 

(or M. N ) whose matrix valued correlation kernel is given by 

E(f (x, s) ® f*(y, s')) = K(x, y, s - s') . (2) 
We will usually assume that K(x, y, t) vanishes for large t, say \t\ > t > 0. 



1.2 Examples 

1.2.1 Schrodinger equation 

Let X be a smooth Riemannian manifold with Laplace-Beltrami operator A. Let 
us give o:I^la non negative function, V a smooth real valued function on 
X, h a non negative constant, and K = —h 2 A + V(x), and take: 

h 

— {u t + a(x)u) + Ku = g . 
i 

It is a particular case of Equation where H = jK + a(x) and / = ^g. Let 
us note for future use that, if H — > 0, the principal symbol of our equation is 
oo + ||£|| 2 + V(x) and a(x) is only a sub-principal term entering into the transport 
equation, but not in the classical dynamics. 
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1.2.2 Wave equations 

Let us start with 

Utt + 2au t — Au = f, a > (3) 

(with A the Laplace-Beltrami operator of a Riemannian metric on X) which 
corresponds to Equation (JJJ with 



. u, J ' " V / 

and 

-Id 

-A 2a 

1.2.3 Pseudo-differential equations 

We can assume that the dynamics is generated by a \PDO (see Section EJ). Our 
equation looks then like: 

-u t + H £ u = f 

% 

with 

£ £ = Op £ (if + £#i) • 

This allows to include 

• An effective surface Hamiltonian associated to stratified media (included in 
the Ho term) j^J. They are usually \I/DO's with a non trivial dispersion 
relation 

• Frequency dependent damping included in H\\ this is usually the case for 
seismic waves. 

1.3 The correlation 

Definition 1 Let us define, for t > 0, Q(t) := exp(-tH) and the propagator P 
by the formula: 

(fi(t)v)(x)= / P(t,x,y)v(y)\dy\ . 
Jx 

The propagator P satisfies 

/ P(t,x,y)P(s,y,z)\dy\ = P(t + s,x,z) 
Jx 

which comes from: Q(t + s) — Q(t) o Q(s). The causal solution of Equation (TJQ) 
is then given by 

poo i* 

u(x,t)= ds P(s,x,y)f(t - s,y)\dy\ . 
Jo Jx 



We define 

1 f T 

C A>B (r) := lim — / u(A,t)®u(B,t-T)*dt . 

Ergodicity allows to replace time average by ensemble average and we get by 
a simple calculation: 

Theorem 1 If P is defined as in Definition^ and K by Equation we have, 
for t > : 

poo rs-\-T p 

C a ,b( t )= ds da \dx\\dy\P(s + T-(T,A,x)K(x,y,cr)P*(s,B,y) 

JO J-oo JXxX 

(4) 

and Ca,b{~ t ) = C b ,a( t Y- 

We get, for r > t > the formula 2 : 

C AB (r) = [n(r)Il}(A,B) (5) 

with 

n = / n(s)cn*(s)ds (6) 

Jo 

and 



C= f Q(-t)K(t)dt . 

J\t\<tn 



\t\<t 

If we assume that K(x,y,o~) = L(x,y)S(a) we get the simpler formula 

C = L . 



2 Exact formulae for white noises 
2.1 Vector valued white noise 

Let us assume that L(x,y) = 5{x — y)ld meaning that f is vector valued white 
noise on I x R. We get, for r > 0, 

poo 

n= / VL(s)VL*(s)ds 
Jo 

and, assuming H = %A + k, with A self-adjoint, the following simple formula: 

C a , b (t) = ±-P (t,A,B) , 

which is an exact relation between the correlation and the propagator P$ of the 
wave equation without attenuation (k = 0). 

2 If R is an operator we will denote by y) its Schwartz kernel; L is the integral operator 

whose kernel is L(x,y) 
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2.2 Twisted white noise 

This section was motivated by a question of Philippe Roux. 

Definition 2 A twisted white noise is a random field given by f = L w with w 
a white noise as defined in Section \2. 1\ and L G £(C , C ). Its correlation is 
5(x — y)S(s — s')K with K = L Lq. 

We have then 

C iiB (r) = [fi(r)n](AB) 

with 

/•oo 

n = / tt(s)K n*(s)ds . 



JO 

In the particular case of the scalar wave equation with a constant damping a 
and the dynamics described in Section ESI by Equation 0, we get using a gauge 
transform where 

n = T^O.diag + R 

where -Ko,diag is the diagonal part of Kq and R is going to vanish in the high 
frequency limit. 

2.3 Wave equations 

Let us take the case of a scalar wave equation with constant damping; the closest 
results were derived in [221 El- 

We will consider the wave equation © with 

• a > is a constant damping coefficient 

• A a Riemannian laplacian in some Riemannian manifold X, possibly with 
boundary: 

A = g ij (x)dij + bi(x)di 

which is self-adjoint with respect to \dx\ and appropriate boundary condi- 
tions; in fact we could replace the Laplacian by any self-adjoint operator 
on X\ 

• f = f(x,t) the source of the noise which will be assumed to be a scalar 
white noise (homogeneous diffuse field): 

E(f(x,s)f(y,s'))=5(s-s')5(x-y) 
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Let us compute the "causal solution" , i.e. the solution given by u = Gf with 
G linear and satisfying u(.,t) = if f{-,s) vanishes for s <t. 
We introduce the vector 



u 



u 



dtu + au 



which satisfies: 



with 



d f \i + au + Hu 



(7) 



H 



Id 

A + a 2 



In order to get a readable expression, it is convenient to introduce P 
y/—(A + a 2 ). We get then easily: 



u(x, t) 



poo r 

/ ds / e~ as 


sin sP 


'o J X 


P 



{x,y)f(y,t- s)\dy\ 



where sin sP is defined from the spectral decomposition of P (any choice of the 
square root gives the same result for sm ^ p ) . The meaning of the brackets is 
"Schwartz kernel of" . We define 



G a {t,x,y) = Y(t) 



-at 



sintP 
P 



with Y the Heaviside function. We will call G a the (causal) Green function. We 
can rewrite: 

u(x,t)= ds G a (t - s,x,y)f(y,s)\dy\ . 



x 



Let us assume now that / is an homogeneous white noise and compute the 
correlation Ca,b(j) := lim^oo -f f u(A,t)u(B,t — r)dt. 
We get quite easily, using 



or: 



sin a sin/3 = -(cos(a — (3) — cos(a + (3)) : 



Ca. b (t) 



-a\r\ 



cos tP a cos rP — P sin It \P 



AaP 2 



Ca,b(t) 
Some comments: 



4(P 2 + a 2 ) 



AP 2 (P 2 + A 2 ) 
cos tP sin \t\P 



(A,B) 



P 



(A,B) 



(8) 
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Low frequency filtering: from Equation JSJ), we see for each eigenmode 
Auj = uJjUj a prefactor, which in the limit of large r is ~ e^^ 3 "^ for 
Uj = a 2 — r 2 < a 2 . This acts as low frequency filter as observed in |13j . 

High frequency regime: in the high frequency regime, we get a simplified 
expression: 

CaAt} ~ 6 4aP 2 
of which, taking the derivative w.r. to r, we get 

„ , s e~ a|r| / sin tP 



4a V p 

recovering the expression derived by many authors: 

d T C A , B (r) « ^ (-G(A, P, r) + G(A, P, -r)) (9) 
4a 

If the attenuation is small, we get 

1 costQ 



Ca,b(t) 



4a Q 2 



with Q = \/—A, and hence the relation (jHJ). 

3 Time reversal symmetry 

Definition 3 • A dispersion relation £, cj) = zs said to fre time reversal 
symmetric (TRS) if it is invariant by a : (x,£) — > (x, — £)■ 

• t4 linear wave equation (with no attenuation! ) is said to be time reversal 
symmetric (TRS) if for any solution u(x,t) the field u(x,—t) is also a 
solution. 

Example 3.1 • Schrddinger equations without magnetic fields 

• Acoustic and elastic wave equations. 

Lemma 1 If D is TRS and j(t) = (x(t),£(t),uj ) is a solution of Hamilton's 
equations, a(j(—t)) = (x(—t),—£(—t),Uo) too. 

We have the following: 
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Proposition 1 The correlation satifies the following general identity: 

C a>b (t) = C* B)A (-t) 

and, in case of a white noise and a time reversible wave dynamics modified by a 
constant attenuation (as in Section Wl\) : 

Ca,b(-t) = C AB (r) ■ (10) 

Approximations of Equation (jlOj) turn out to be important in applications to 
clocks synchronisation. 

4 Random fields and pseudo-differential opera- 
tors 

4.1 Goal 

Our aim in this section is to build quite general random fields with correlation 
distances given by a small parameter e. It seems to be natural for that purpose to 
use e— pseudo-differential operators. We will see how to compute the generalized 
power spectrum using Wigner measures. 

4.2 White noises 

Let (7i, (.|.)) be an Hilbert space. There exists a canonical Gaussian random 
field on it, called the white noise and denoted by (or simply w if there is no 
possible confusion). This random field is defined by the properties that: 

• For all eeH, 

E((w|e)) = 

• For all e,f G H, 

E((w|e)(w|/)) = (e1/) 

Unfortunately, w is not a random vector in 7i unless dim 7i < oo 3 , but only 
a random Schwartz distribution. 

We have nevertheless the following usefull proposition: 

Proposition 2 If A is an Hilbert- Schmidt operator on 7i, the random field Aw 
is almost surely in Tt. 

3 If w were a vector in 7i, we would have w = XX w l e i) e i f° r an y orthonormal basis (€j) 
and we see that 

E(||w|| 2 ) = ^E((w|e-)) 2 =dimH. 
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Proof. - 

E({Aw\Aw)) = E((A*Aw|w)) = Trace(AM) which is finite, by 
definition, exactly for Hilbert-Schmidt operators. 

□ 

4.3 Examples 

Example 4.1 Stationary noise on the real line: let us take a random field on the 
real line which is given by the convolution product of the scalar white noise w with 
a fixed smooth compactly supported function F: f — F-kw. Then f is stationary: 
it means that the correlation kernel K(t,t') = E(/(i)/(i')) is a fonction oft — t' . 

On the level of Fourier transforms f = Fw, E(/(a;)/(a; / )) = \F\ 2 (lj)5(lj — uo 1 ) and 
the positive function \F\ 2 (uj) is usually called the power spectrum of the stationary 
noise. 

Example 4.2 If X is a d-dimensional bounded domain. Let us denote the Sobolev 
spaces on X by H S (X). If P : L 2 (X) -> H S (X) with s > d/2, Pw is in L 2 (X). 

Example 4.3 Brownian motions: if X = R ; w is the derivative of the Brownian 
motion: if b{t) = J^w(s)ds, b : [0,+oo[— > R is the Brownian motion which is in 
L 2 ([0,T]) for all finite T. 

Example 4.4 If X is a smooth compact manifold or domain and P is smoothing, 
meaning that P is given by an integral smooth kernel 

Pf(x)= [ [P](x,y)f(y)\dy\ , 
Jx 

F = Pw is a random smooth function. Its correlation kernel 

C{x,y):=E{F(x)F(y)) 

is given by: 

[PP*](x,y) = [ [P](x,z)[P](y,z)\dz\ . 
Jx 

Example 4.5 Random vector fields: let us consider 7i = L 2 (X, ~R N ). For ex- 
ample, in the case of elasticity, X is a 3D domain and N = 3. The field here are 
just fields of infinitesimal deformations (a vector field) . 

4.4 Modelling the noise using pseudo-differential opera- 
tors 

The main goal of the present section is to build natural random fields which are 
non homogeneous with small distances of correlation of the order of e — ► 0. The 
noise is non homogeneous in X, but could also be non isotropic w.r. to directions. 
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4.4.1 Noises from pseudo-differential operators 



It is therefore natural to take for noise on a manifold Z the image of an homo- 
geneous white noise by a pseudo-differential operator N of smooth compactly 
supported symbol n(z,Q. The correlation C(z,z') will then be given as the 
Schwartz kernel of NN* which is a ^>DO of principal symbol \n\ 2 . We have 



while H 2 (z, C) is the "power spectrum" of the noise at the point z. 

This construction gives smooth random fields which can be localized in some 
very small domains of the manifold Z, which are non isotropic and which have 
small distance of correlations. Moreover it will allow to use technics of microlocal 
analysis with the small parameter given by e. 

4.4.2 Power spectrum and Wigner measures 

Definition 4 If f = (f £ ) is a suitable family of functions on Z , the Wigner 
measures W e f of f are the signed measures on the phase space T*Z defined by 



The measures dWf are the phase space densities of energy of the functions f £ . 
We now define: 

Definition 5 The power spectrum of the random field f = (f £ ) is the phase 
space density PI defined by: 





P £ f = E(W £ f ) 



the power spectrum of a random field is its average Wigner measure. 
Proposition 3 The power spectrum P of f £ = Op £ (n)w, satisfies: 



P e f ~ (2TXEY d \n\\x^)\dxdi\ . 



Proof. - 



Let us put N = Op £ (ri), we have 



(Op £ (a)Nw\Nw) = (N*Op £ (a)Nw\w) 



and E((Aw\w)) = trace(A). We get 




which can be evaluated using the \l/DO calculus as 




□ 
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4.4.3 Space-time noises 

If Z = X x R is the space-time, we will take our noise as before / = Lw; we will 
assume the noise homogeneous in time, the symbol I of L is assumed to be given 
by l(x,Z,w). 

In this case, the correlation is given by: 

K(x,y;t) = [LL*](x,y;0,t) (11) 

which is the Schwartz kernel of a ^>DO of principal symbol 11* (x, £; oS). 

4.5 Equipartition and polarizations 

If f e : Z — > R 9 is a vector valued family of functions, the Wigner measure is 
matrix valued: if a : T*X — > Sym(R 9 ), we define 

j adWt := (Op £ (a)f £ |f £ ) 

If we have, for each point (z, C) of T*Z a splitting R 9 = ®Ej(z, () with projectors 
Pj, there exists canonical measures fij defined by 

J adfij = Trace (aPj)\dzd(\ . 

If the Ej's are defined by the polarizations of an Hamiltonian H, the micro- 
canonical Liouville measures are given by 

/adL-E := / ad/Xj . 
j J\j<E 

A (random) state of energy E of H is said to equipartited if its (average) 
Wigner measure converges to the microcanonical measure. If P% is the spectral 
projector of H over the modes of energies less than E, the random state Pew is 
equipartited. 



5 High frequency limit of the correlation: Schrodinger 
equations 

The main result is easier to derive in the case of a scalar field gouverned by a wave 
equation which gives the first order derivative of the field: it is a generalization 
of the Schrodinger equation. 
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5.1 Assumptions 

Let us start with the semi-classical Schrodinger like equation 

e - e 
-u t + H e u = -f 

i % 

where 

• H e = ~H is an e— pseudo-differential operator: 

# e :=Op £ (# + £#i) 

with 

— The principal symbol Hq(x,£) : T*X — > K. which gives the classical 
("rays") dynamics: 

dxj = dHo d^i = 1 < . < d 

dt d£j ' dt dxj ' ~~ ~~ 

We will denote by 4> t the flow of the previous vector field. 

— For technical reasons (see Lemma OJ), we assume that Hq is elliptic at 
infinity: lim^oo H (x, £) = +00. We define Ho = Op e (H ) and the 
unitary group U(t) = exp(—itH /e). We define also H\ = Op e (ifi). 

— The sub-principal symbol Hi(x,£) admits some positivity property 
which controles the attenuation: there exists k > 0, such that 

Sifi < —k . 

• The random field / is given by / = Op e (l(x,£,u))w with w the white 
noise on X x K and with / smooth, compactly supported w.r. to (x,£) 
and whose Fourier transform w.r. to uj is compactly supported. The power 
spectrum of / is (27r£r)~^ +1 ^| 2 (x, £, uj). 

The previous assumptions will be used everywhere inside Sectional 

5.2 Subprincipal symbols and attenuation 

Lemma 2 Under the assumptions of Section \5.1\ we have: 

a) There exists c> so that, for \t\ < c\ log(e:)|, fl(t) = U(t)Y(t) with Y(t) a 
of principal symbol exp(— i J Q Hi(<f) s (x,£))ds). 

b) for allt > 0, the estimates 

\\Y(t)\\=0(e- k,i ) 
with \\.\\ the operator norm in L 2 (X) and < k! < k. 
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Proof. - 

a) We start with fi(t) = U(t)Y(t) and hence 

Y'(t) + iUi-QHxUtyYQ) = . 

Using Egorov's Theorem with logarithmic times, as in Section E3J we 
get 

Y'(t)+iH 1 (t)Y(t) = 

with Hi(t) a \&DO of principal symbol iHi(4> t (x,£)). It is then enough 
to start with a formal expansion in e of the symbol of Y(t) and to 
work by induction on the powers of e. 
b) We have 

j t (v(t)\v(t))^2^(v(t)\-iH 1 (t)v(t)) 

and we use Garding inequality (see [7j): if a > 0, Op £ (a) > — C for 
any C > and £ small enough. 

□ 

5.3 Some lemmas 

Lemma 3 If A = Op e (a) with a compactly supported, the operator B = exp(itH)A 
is a ^/DO of principal symbol b = exp(itH )a. 

Proof. - 

Let us choose a function \ £ C^°(R, R) so that x(i?o) is equal to 1 
in some neighbourhood of the support of a. We have 

B = exp(itH ) (x(H ) + (1 - x(#o))) C 

with C = Y(—et)A is a compactly supported \l/DO of principal 
symbol a. The previous expression of B splits into 2 terms I? = I +11. 
The first one rewrites I = <&(Hq)C with $ = e lt '\ £ which is a 
^DO of principal symbol e °a thanks to the functional calculus 
of elliptic self-adjoint $DO's (see [7j, Chap. 8). The second one is 
smoothing. 

□ 

The following Lemma follows directly from the definitions: 

Lemma 4 If K(x,y,t) is the correlation kernel of f , then K(ai), the operator 
whose kernel is K(., ., eai), is a ^>DO which vanishes for \a\\ > C and whose 
principal symbol is {2tis)~ 1 j e lr71UJ \l\ 2 (x, £, u)du. 
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From Equation we get, for r > and e < C/t: 

poo 

C A:B (t) = [fi(r) / dsfi(s)£fi*(s)](A,fl) (12) 

with 

£ = e Q(—eai)K(ai)dai . 

J\(Tl\<C 

Lemma 5 £ is a ^DO of principal symbol \l\ 2 (x, £, —H (x,£)). 
Proof. - 

The result follows from Lemma El and the value of the symbol of 
K(ai) given in Lemma EJ by integrating w.r. to a. 

□ 



5.4 Applying Egorov Theorem 

We can apply Egorov Theorem: 

Theorem 2 (Egorov's Theorem) If A = Op £ (a), U(-t)AU(t) = A t where A t 
is a ^>DO of principal symbol a o <f> t with <p t is the Hamiltonian flow of H . 

We will need a large time estimation in the Egorov's Theorem. Such estimates 
are provided in the nice paper £Q: Egorov's Theorem still works under suitable 
hypothesis on Hq for time bounded by c| loge| where c is related to the Liapounov 
exponent of the classical dynamics. Such time is called Eherenfest time and will 
be denoted by T Ehrenfest . 
We get the main result: 

Theorem 3 With the assumptions of Section \5.1l the correlation is given, for 
t > 0, by 

C a>b (t) = P(t)oIL)](A,B) 
where LT = Op e (7r) + R with: 

ir(x,£) = f exp (2 r^H^ix^ds) \l\ 2 0, ~ H o(^ 0) dt • 

J~c\\oge\ \ Jt J 

and R the remainder term is u O{e a )". More precisely, let us consider C 'a,b{ t ) 
the Schwartz kernel of an operator C(r). This operator is Hilbert- Schmidt 4 with 
an Hilbert- Schmidt norm of the order of e~ d l 2 . We have 

||£||h-s = 0{e a ~ d ' 2 ) . 

4 An Hilbert-Schmidt operator A is an operator whose Schwartz kernel [A](x, y) is in L 2 (X x 
X) and the Hilbert-Schmidt norm ||A||h-s of A is the L 2 norm of [A]. 
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Proof. - 

We start from 

/■oo 

11= / tt(s)£n*(s)ds (13) 
Jo 

as given by Equations (HJ) and (fT^j) . We have = ^( s )^( s ) w ^ n 
y(s) is a \l/DO of principal symbol exp(— i J" Hi{(f) u {x^))du) and 

Q(s)£fi*(s) = U{-s) 

to which we want to apply Egorov's Theorem. There is a techni- 
cal problem due to the fact that Egorov's Theorem is only valid 
untill Ehrenfest times. We split the integral (|13j) into two parts 
II = L Ehrenfost + f°° . The first term is estimated using Egorov's 

J U J -L Ehrenfest 

Theorem for large times. An upper bound for second part follows 
from the decay estimate of gi ven m Lemma El and the esti- 

mates ||Op e (a)|| H -s = 0(e~ d / 2 ). 

□ 

Assuming still r > 0, we see that the correlation Ca,b(j) is close to the kernel 
of a Fourier integral operator associated to the canonical transformation <p T . It is 
given as a sum over all classical trajectories 7 from B to A in time r of Cauchy 
data (B,£b) with H (B,£b) = f° r which the backward trajectories crosses 
the support of the power spectrum //*(., ., to). If 7 is such a trajectory and B and 
A are non conjugated along it, this contribution is given by the well known Van 
Vleck formula 5 multiplied by tt(B,^b)- 

Corollary 1 Let K be the support of —Hq(x,^)) and the smallest 

closed set of T*X invariant by the Hamiltonian flow of Hq and containing K . 
The Hamiltonian Hq restricted to can be recovered from the knowledge of 
C{t) for0< |t| < r . 

In particular, if there exists (x,(,) with H (x,£,) = E and l(x, £, —E) 7^ and 
if 4> t is ergodic on H^ l {E), then we can recover the flow <p t on H^^E). 

5.5 Remarks 

We would like to extend the previous approach to the general case, i.e. to iV > 1. 
There are 2 difficulties to overcome: 

• One has to extend Lemma El to the case of systems 

• Egorov Theorem is no more true, but remains true on average as in Lemma 

H 

5 The Van Vleck formula expresses the propagator P(t, A, B) as a sum of p 7 = 
(2Trie)~ d ^ 2 a 1 (s)exp(iS( , y)/e) with a 7 (e) a formal power series in e with a first term explicitely 
computable 
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6 High frequency limit of the correlation: wave 
equations 

We want to derive results similar to those of Theorem El in the case of the scalar 
wave equation (N = 2) given as follows: 

u tt + 2au t - Au = f 

where a > is constant and —A is the Laplace-Beltrami operator on a smooth 
complete Riemannian manifold X. 

We will assume that the source noise / is given by / = Lw with w a white 
noise onlxR and L = Op e (Z) with I = smooth, with compact support, 

and independant of to. We will moreover assume time reversal symmetry of the 
noise, namely l(x, — £) = l(x, £); it implies that the kernel [L](x,y) is real valued. 



6.1 Direct derivation 



Let us introduce P = a/— A — a 2 , the causal solution of Equation Q is given, as 
already used in Section 12.31 by 



u(x, t) 

with the Green function 



ds I G(s, x,y)u{t — s,y)\dy\ 



x 



G(t,x,y) = Y(t)e 



-at 



sin tP 



By a direct calculation and denoting by K = P l LL k P 1 with L = Op e (Z), we 
get, for r > 0, 



Ca,b{t) 



-e 
2 



3? e 



-irP 



e 2as e isP 



K (e~ is 



') ds 



(A,B). (14) 



This integral splits into 2 parts, the first one can be asymptotically computed 
using Egorov's Theorem as in the proof of Theorem El while the second is smaller 
by the 

Lemma 6 If A is a compactly supported ^>DO of order and if 

r-0 



J :-- 



e 2as e isP Ae isP ds 



we have 

with some non negative 7. 



|J|| H -s = 0(e 
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This is proved by cuting the integral into 2 pieces as in the proof of Theorem EJ 
using Egorov's Theorem and integrating by part. 
The final result is: 

Theorem 4 With the previous assumptions, we have, for r > 0: 



and \\R\\ns = 0{e^ +2 - d ' 2 ) with 7 > 0. 

Remark 1 The pref actor e 2 is just here because we have to pass from the ^DO's 
without small parameter P to an e—^DO: P = £ -1 Op e (||£||) + l.o.t.. 

6.2 Using the general formalism 

We can again start from Equation (@). We are reduced to calculate the correlation 
between u(A,t) and u(B,t) which is given by 



We should first put the wave equation as a first order semi-clasical equation. 
We put 



— e~ ar [cos rP o n] (A, B) + R 



with n a tyDO of principal symbol 




C 1 A ] B (r) = Mr)U)Y\A,B) . 




u 

e{u t + au) 




and 




We get, by definig P = a/— (A + a 2 ) 




Moreover, if / = Op e (m)w, 



f 




with 




We see that the computation is less easy • • • 



18 



7 What's left? 



There are still several problem to discuss: 

• The case of vector valued wave equations with several polarizations. 

• The precise study of autocorrelations: as mentionned to me by U. Smilanski, 
the autocorrelations can be usefull in order to learn more about the source 
of the noise f . 

• The case of surface waves: effective Hamiltonains and the associated inverse 
spectral problems will be discussed in [5]. 

• The case of the source noise located on the boundary: 



will be discussed in jO]. 

8 Random scattered waves 

In this last independant section, we will revisit what was maybe the starting 
point of this story by Keiiti Aki in the fifties: he wanted to measure the speed of 
propagation of seismic plane waves by averaging over the incidence directions. It 
turns out that we get nice formulae even for non homogoeneous media. 

8.1 Introduction 

Let us consider the propagation of waves outside a compact domain D in the 
Euclidian space M d . Let us put fl = M. d \D. We can assume for example Neumann 
boundary conditions. We will denote by Aq the previous self-adjoint operator. 
So our stationary wave equation is the Helmoltz equation Anf + k 2 f = with the 
boundary conditions. We consider a bounded intervall / = [E_,E + ] c]0, +oo[ 
and the Hilbert subspace Tij of L 2 (Q) which is the image of the spectral projector 
Pi of our Laplace operator Aq. 

Let us compute the integral kernel 11/ (x, y) of Pi defined by: 



u tt + au t — An = 
UdX = f 




into 2 different ways: 



1. From general spectral theory 



2. From scattering theory. 
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Taking the derivatives of Ui(x,y) w.r. to E + , we get a simple general and ex- 
act relation between the correlation of scattered waves and the Green's function 
confirming the calculations from ^1] in the case where D is a disk. 



8.2 Hj(x,y) from spectral theory 

Using the resolvent kernel (Green's function) G(k,x,y) = [(k 2 + Aq)~ 1 ](x,|/) for 
Qk > and the Stone formula, we have: 

n 7 (x, y) = ^(^JJ G{k + iO, x, y)kdk\ 

Taking the derivative w.r. to k + of U[E^,k 2 ]( x , v), we get 
d 2k 

M^[E-^] (x, y) = —%G(k + iO, x, y)) . (15) 

8.3 Short review of scattering theory 

They are many references for scattering theory: for example [T2*] . 
Let us define the plane waves 

e (x,k) = e i<k ' :l ' > . 

We are looking for solutions 

e(x, k) = eo(x, k) + e s (x, k) 

of the Helmoltz equation in Q where e s , the scattered wave satisfies the so-called 
Sommerfeld radiation condition: 

^ k) = -^j- 2 ( e°°(o, k) + O(o)l ,x^oo. 



The complex function e°°(x, k) is usually called the scattering amplitude. 

It is known that the previous problem admits an unique solution. In more 
physical terms, e(x, k) is the wave generated by the full scattering process from 
the plane wave eo{x,k). Moreover we have a generalized Fourier transform: 



f(x) = (2ny d [ /(k)e(x,k)|dk| 

JR d 



with 



/(k)= / e(y,k)f(y)\dy\ 



From the previous generalized Fourier transform, we can get the kernel of any 
function $(— Aq) as follows: 

M-A n )](x,y) = (2ny d [ $(A; 2 )e(x, k)i(^k)|dk| . (16) 
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8.4 Uj(x,y) from scattering theory 

Using Equation ([TH| with $ = 1/ the characteristic functions of some bounded 
intervall /, we get: 

n 7 (x, y ) = ( 27 r)- d / e(x, k)^k) \dk\ . 

Je_<u. 2 <e + 

Using polar coordinates and defining \da\ as the usual measure on the unit (d — 
1)— dimensional sphere, we get: 

n 7 (x,y) = (27r)~ d f [ e(x,k)^k)k d - l dk\da\ . 

We will denote by <r d _i the total volume of the unit sphere in R d : a = 2, o x = 
2n, cr 2 = 4-7T, ■ ■ • . 

Taking the same derivative as before, we get: 



This integral can be interpreted, using the correlation Cf? att (a;, y) of random 
scattered waves of energy E defined by 

CT U (x,y) = — f e(x,k)^k)\da\, 

&d-l Jk 2 =E 

as 

^ {E -H^y)-^^cr\^v) ■ (17) 

8.5 Correlation of scattered plane waves and Green's func- 
tion: the scalar case 

From Equations (fTKJ) and (fTTjl . we get: 

^±lCT U (x,y) = ^(G(k + iO,x,y)) . 

Hence 

For later use, we put 

nd+l d-l 

- • (18) 
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8.6 The case of elastic waves 

We will consider the vectorial stationary elastic wave equation in the domain Q: 

Hu - uo 2 u = 0, 
with symmetric boundary conditions, where 

Hu = —a Au — b grad divu . 
where a and b are constant: 

H , A + /i 
a = — , o = 

P P 

with A, n the Lame's coefficients and p the density of the medium. 

• The case n = R d 

We want to derive the spectral decomposition of H from the Fourier inver- 
sion formula. Let us choose, for k ^ 0, by k, k 1; • ■ ■ , kd-i an orthonormal 
basis of M. d with k = | such that these vectors depends in a measurable 
way of k. Let us introduce Pp = kk the orthogonal projector onto k 
and Pg = Yl'jZi kjk* so that Pp + Ps = Id. Those projectors correspond 
respectively to the polarizations of P— and S— waves. 

We have 

ILfay) = (270-^^-^ ((a + b)~ d / 2 J k2=u)2/(a+b)2 e i ^P*da+ 
using the plane waves 

e£(x,k) =e aix k 

and 

e s,ji x i k) = e* x kj 

we get the formula: 

Uj(x, y) = (27r)- d l 2eI u d ^du ((a + b)^ J k2=uj2/(a+b)2 e?(s, k)(e°(y, k))*da+ 
«" d/2 E -I! 1 L^ /a > e° s>j (x, k)(e° SJ (y, k)) W) . 

• Scattered plane waves 

There exists scattered plane waves 

ep(x, k) = e P (x, k) + e P (x, k) 

e Stj (x, k) = e^(z, k) + e s S j (x, k) 

satisfying the Sommerfeld condition and from which we can deduce the 
spectral decomposition of H. 
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• Correlations of scattered plane waves and Green's function 

Following the same path as for scalar waves, we get an identity which holds 
now for the full Green's tensor QG(w + iO, x, y): 

3G(w + iO, x, y) = -y d (u) ({a + b)~ d / 2 J k2=uj2 /(a+b)2 e P (x, k)(e P (y, k))*da + 
with 7d(u;) defined by Equation fT8|) . 

This formula expresses the fact that the correlation of scattered plane waves 
randomized with the appropriate weights is proportional to the Green's 
tensor. 

9 Appendix: A short review about pseudo-differential 
operators 

We will define pseudo-differential operator ($DO's) on M. d . ^DO's on manifold 
are defined locally by the same formulae. More details can be foud in |HI71IHIIT^]. 

Definition 6 • The space is the space of smooth functions p : T*M. d — > C 
which satisfies 

Va,/3, \D«Dlp(x,0\ < Cajil + |e|) fc " l/31 • 

• A symbol of order m and degree I is a family of functions 

p £ : T*M. d -> C 
which admits an asymptotic expansion 

oo 

p E ^Y, £j+m p^o 

3=0 

with pj G . We will denote this space by S m j. 

Definition 7 An e-pseudo-differential operator P (a ^DO ) of order m and de- 
gree I on lR d is given locally by the kernel 

[P](z,z') = (27T6)- rf l/ (Z - Z ' l0/E Pe (^>C) Kl 

where p e {z,Q, the so-called (total) symbol of P, is in S m j. 
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We will denote P = Op £ (p £ ). 
The kernel of P is then given by: 

with p the partial Fourier transform of p e (z, () w.r. to (. Very often, one is only 
able to compute the symbol po which is called the principal symbol of P. 

The most basic fact about MO's is the fact they can be composed: if P = 
Op e (p) and Q = Op E (q), we have PQ = Op e (pq + 0(e)). 
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